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1. Introduction 

Formation of the confining string between external heavy quark and antiquark naturally ex- 
plains the linearly rising heavy quark potential as well as specific stringy corrections present in 
it This phenomenon is confirmed in numerous lattice studies (see, e.g., Refs. [^), which how- 
ever are always limited to some finite range of UV cutoff (lattice spacing a). However, the internal 
structure of the confining string in terms of various gluonic observables is less known, although it is 
generically expected and confirmed on the lattice |Q, |3|] that within the flux tube the energy density 
is enhanced, while the action density, the magnitude of topological fluctuations and other quantities 
are suppressed with respect to the corresponding vacuum expectations. Note that the local action 
density difference 

A.= (.)o-(.),^>0, s = TrF^^ (1.1) 

is distinguished since its integral matches the heavy quark potential [^|. In the Umit of large qq 
distances R the corresponding relation for SU(2) gauge group reads 

d^rAs = 2aP„ ^/ = -^-|^> (1-2) 

where the integral is over the flux tube transverse cross section and a is the string tension. Note that 
the rigorous sum rule (|r|) says hardly anything about actual action density profile. Nevertheless, 
one expects rather generically that its functional form is to be Gaussian [§] 

As{r,R) = C{R)cxp{-r^/d\R)}, 5\R) = X ln[/?//?o] , (1-3) 

where we quoted common to all effective string theories prediction for the squared width 5^ of the 
flux tube, is model-dependent numerical constant and Rq is UV cutoff of the effective theory. True 
that the string action density profile attracted much attention in the past [^. The summary of these 
investigations is that both the Gaussian shape and the string widening, Eq. ( |1.3[ ), are qualitatively 
confirmed for large enough quark-antiquark separations R > 0.7 fm. Although quantitatively there 
are points worth to be improved, the above qualitative picture will hardly change. Thus the question 
is why we decided to consider the flux tube profile again and what are the subtle points we hope to 
shed some light on? 

Briefly, the reason reduces to the recent observation that non-perturbative fluctuations, com- 
monly considered in relation to confinement and known as Abelian monopoles and center vortices, 
have a non-trivial ultraviolet structure |^ (see |^ ||] for comprehensive discussions). Namely, the 
excess of non-Abelian action associated with both the monopoles and vortices is UV divergent, 
while their number density perfectly scales towards the continuum limit. On the other hand, in 
Yang-Mills theory all UV divergences must be interpretable within the field theory itself. Con- 
fronting the field theory with the lattice data one concludes that vacuum defects (monopoles and 
vortices) are, in a sense, dual to high orders of perturbation theory (see [Q] for details). For the 
present discussion it is crucial that for the action density the impact of the vacuum defects (or, 
alternatively, of the ambiguities of the perturbation theory at high orders) is to be parametrized by 
an additional quadratic term in the conventional OPE-based expansion 

(^)o = • I L < + (a • Aqcd? + C4 (a • Aqcd)^ I . (1.4) 



2 



Width of the Confining String 



F.V. Gubarev 



The necessity of the quadratic correction is quite evident in the context of vacuum defects. How- 



ever, from the field-theoretical viewpoint, which is actually implied by ( |1.4| ), a remark have to be 



added. Namely, Eq. ( |1 .4| ) explicitly refers to a particular N-loop order of perturbative expansion. 



the divergent reminder of which is parametrized by the second term. In this respect ( |1.4| ) does not 
follow from the theory itself, rather it reflects the inability to handle high orders of a^-expansion 
and the successful phenomenology of quadratic power corrections (see Ref. ^] and references 



therein). In particular, by no means Eq. ( |1.4| ) implies OPE violation. Consequently, in the field 
theory settings the quadratic correction cannot be discussed out of context of the perturbation the- 
ory. This is reflected in the explicit A/^-dependence of the coefficient : with increasing number 
of perturbative loops explicitly accounted for it drops rapidly being (numerically) compatible with 
zero in the empirical limit N —>■ °° [p^. Nevertheless, for moderate N, when the quadratic term 



could be detected, it turns out to be parametrically small QllQ 



1-2 



2 < 

QCD 



[50MeV]^ 



(1.5) 



which is only an order of magnitude estimation (concrete numbers vary depending on the details 
of the measurements). 

Given the above uncertain status of the quadratic divergence it is legitimate to ask whether the 
issue is at all relevant. The fair answer seems to be affirmative because we aware of the signifi- 
cance of vacuum defects to QCD vacuum structure. The discovery of non-trivial UV structure of 
monopoles/vortices provides a bridge to the field theory Hnking them to the question of power di- 



vergences. Hence it is worth to consider Eq. ( |1.4[ ) in its generality and to consider the most general 
power dependence of {s)g upon the UV/IR scales 



2 + r(°) 



(1.6) 



with some unknown coefficients. True that Eq. ( |1.6[ ) by itself makes no much sense because the 
vacuum action density is not observable anyway. However, it is crucial that once the quadratic term 
does not vanish identically no symmetry could ever prevent its dependence upon the external soft 
fields. Hence the corresponding VEV {s)^^^ is to be parametrized analogously with a priori distinct 



coefficients 0^2^ • Then the difference ( |1.1| ) inside the string is given by 



^s = C2A; 



QCD 



QCD-, 



(1.7) 



while Eq. (L5) provides an upper bound on C2Ag(^^. It is worth mentioning that in terms of vac- 
uum defects Eq. (L7) with C2 7^ follows directly from the known distribution of monopoles [12] 
and vortices [ pj| ] in the vicinity of the flux tube. Combining (|]^, ( |L7| ) one concludes that non- 
vanishing C2 coefficient manifests itself in peculiar linear shrinkage of the flux tube width with 



simultaneous quadratic divergence of the height of the Gaussian profile (1.3) 



QCD: 



C{R) 



A2 

^^QCD " 



(1.8) 



in the limit a ^ (insisting on Gaussian profile seems to be safe at any finite a). Note that this be- 



havior is only asymptotic, an estimation of the relevant lattice spacings below which ( |1.8[ ) settles in 
could be obtained from (L5) and by noting that the term C4 Agj^^, is to be of order gluon condensate. 
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J3 


a, fm 




Multileve 
Torelons, Niev = 2 


;1: A and N,ipd 
Wilson loops, Niev = 3 


Multichannel, 

^smearing 


2.5046 


0.081 


28'^, [2.3 fm]'^ 


(14), (4000) 






2.5515 


0.072 


28"*, p.Ofm]'* 
32^, [2.3 fm]'^ 


(14), (3500) 
(16), (5000) 


(4,8), (12,300) 


(50,65,80) 


2.6800 


0.047 


36*, [1.7fm]4 
40"*, [l.Qfm]'^ 


(18), (4500) 
(20), (5000) 


(5,10), (12,600) 


(150,180,210) 


2.7600 


0.034 


48^ [1.8fm]'^ 


(24), (5000) 


(6,12), (12,1200) 


(180,210,240) 



Table 1: Ensembles used in our study. Multilevel notations (fourth column) closely follow [Q. APE 
smearing (fifth column) had been performed with £ parameter 0.25. 



Straightforward calculation shows that the divergence (L8) should be manifest for a < 0.03 fm and 



signifies that numerically the status of quadratic divergence is also uncertain. Indeed, the relevant 



scale is rather small parametrically and still leaves slight possibility that the divergence ( |1.7| ) could 
have been missed in the previous studies. The prime purpose of our investigations, which we sum- 
marize in this paper, is to study the scaling properties of the action density distribution within the 
long confining string. For technical reasons we restricted ourselves to SU(2) pure gauge theory 
with Wilson action, summary of simulation parameters is given in Tabled 



Our first measurements 014[ ] confirmed the presence of the quadratic term, Eq. ( |1.7[ ), with 



magnitude being in agreement with (|1.5|). However, we had never been quite satisfied with this 

result for the following reasons: 

1. Ground state separation (ground state overlap). 



String creation/annihilation operators used in [ 14 1 were not optimal, the relevant overlaps with the 



ground state of heavy ^17-pair were of order 80-^90%. This required us to make rather compUcated 
fit to r ^ 00 thus introducing systematic biases the magnitude of which we were unable to fairly 
estimate. This issue is addressed in Section ^. 
2. Limit of large qq separations. 

Ultimately, we're interested in the physics of long confining string, /? » 7?o ~ 0.7 fm, when the qq 
separation R greatly exceeds the string formation length Rq. With Wilson loops (or Polyakov lines) 
it is hardly possible to go far beyond Rq due to the technical limitations, hence the string profile 
is disturbed by singular point-like sources at its end-points. Again some phenomenology had been 
applied in [0] leading to the same question of systematic errors. Ideally, the point-like sources 
have to be removed, which leads us to consider the torelon correlation function (Section ^. 

An overall conclusion of the present study is that it is the overlap problem which is the most 
crucial to get the physics right. In particular, when the proper selection of string ground state is 
ensured and its T = 00 asymptotic structure is directly seen the effect of string shrinkage disappears. 
To our surprise the string profile turned to be extremely sensitive in this respect, contrary to, e.g., 
heavy quark potential in which we never observed scaling violations. Therefore in string profile 
studies it is absolutely necessary to get rid of excited states contamination prior to making any 
conclusions on the flux tube internal structure. 
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2. Torelons 

Technically, the torelon operator is similar to Polyakov line: it is the gauge holonomy around 
the spatial period of toroidal lattice in a particular time slice (see, e.g., Ref. | [l5[ ] and references 
therein). Physically, it is quite different however: torelon creates confining string closed through 
the periodic boundary conditions and stable for topological reasons. The string exists without point- 
like sources, hence the torelon correlation function is best suited to study the flux tube internals. 
Normalized torelon correlator allows usual spectral decomposition --^{t) = c/e^'"'', ^(0) = 1, 
where masses m,- depend upon the spatial lattice extent L and have the same leading asymptotic 
Mi f» gL. 

In order to study the string profile with torelons it is necessary to have precise data for the 
torelon correlation function at large Euclidean times. This is due to the local action density inser- 
tion: for instance, we used clover-like sum of all neighboring to given point plaquettes, hence the 
minimal Euclidean time we could consider is ? = 4 (in fact, it has to be larger for asymptotic expo- 
nential falloff to take place). Thus the use of modern multilevel technique [1^ 17] is mandatory. 



In usual multilevel approaches it is the time axis which is divided into slabs, so that at large Eu- 
clidean times torelon creation/annihilation operators belong to different time slabs. Our experience 
showed, however, that in order to obtain desired precision in string profile, the 2-point correlator 
should be known with by order of magnitude better accuracy than it is available with usual scheme. 
Therefore, we utilized another approach in which the torelon direction is divided among different 
slabs. The ground state enhancement is done as usual in each slab, however, it is evident that our 
scheme does not allow to get optimal overlap with the string ground state: the more slabs are intro- 
duced the worse is overlap. Indeed, in each slab one considers straight smeared spatial holonomies 
with end-points fixed at the boundaries. Upon averaging the string wavefunction is constructed as 
superposition of various paths, which are required, however, to pierce prescribed boundary points. 
Effectively at each boundary one introduces a node of the wavefunction thus reducing the ground 
state overlap. Minimal number of time slabs with the 2-level scheme is two and the corresponding 
overlap is about 50%. Despite of its smallness we hoped to eliminate exited states contamination 
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by attaining the limit of large Euclidean times, which is feasible with multilevel strategy. 

Sample of data points obtained for the torelon correlator is shown on Fig. [T] (left). Analyzes 
reveal that we see a superposition of two states with approximately equal amplitudes. Mass extrac- 
tion is done via two-exponential fit, lower f-limit of which was varied until lighter mass stabilizes. 
Masses extracted this way are in qualitative agreement with the literature [jl^], however, quantita- 
tively they appear to be systematically slightly larger. This is because we are obliged to consider 
2-point correlator, not the wall-wall correlation function, and hence inevitably get an admixture of 
non-zero lowest spatial momentum states. Thus the masses extracted from 2-point function could 
only be used as an approximation to the true spectral gap in the torelon channel. 

Unfortunately, our hope to observe exclusively the ground state profile at largest available 
times failed for small transverse distances (Fig. [T], right). Indeed, for < 0.03 fm^ f-dependence 
of As{t,r^) is apparent, while at larger an approximate f-constancy of the profile seems to settle 
at f > 0.5 fm. Therefore, it is necessary to make t —>-oo extrapolation for which the leading finite-? 
correction looks like single exponent with gap being fixed in term of the spectral gap in the torelon 
channel. However, besides the imprecise knowledge of the gap there is an additional limitation: 
due to imperfect overlap with the ground state single exponent might not be adequate especially 
at small r. In view of these shortcomings various fits were tried out: 1. Single exponent fits 
with gap fixed either from the 2-point function or from the bosonic string theory prediction (they 
are illustrated Fig. |l], right). 2. Fit by constant in f-range where data points come to plateau, it 
was used for qualitative purposes only. 3. Cumulative fit to all available t, points with single 
exponential ansatz and treating the gap as free parameter. To ensure the validity of single exponent 
approximation the points too close to the string axis were excluded, the corresponding cut on r is 
obtained from previous fits. Additionally, we made seemingly well grounded assumption that the 
string profile at f = oo is to be Gaussian. This greatly reduces the dimensionahty of the problem 
and allows much better precision. 

Summary of the results obtained with various fitting procedures is presented on Fig. |2[ left 
panel. Briefly, all the above fits appear to be adequate and consistent with each other except for the 
region < 0.03 fm^ at smallest lattice spacing, where fixed gap fits notably deviate from Gaussian. 
Treating the discrepancies at small r as a consequence of the imperfect overlap and the enforcement 
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Figure 3: Left: As vs. t at (fo,f ) = {(8, 18), (8,22)}, R = 0.93 fm and different (Ji = 2.680, 40"^) with the 
corresponding fits to constant. Right: flux tube profiles at various lattice resolutions. 



of single exponential anzatz, we conclude that regardless of actual t ^ oo extrapolation the string 
width seems to drop with diminishing lattice spacing. In fact, the apparent string shrinkage could 
be identified directly by looking on the data at roughly the same transverse distance r ?a 0.27 fm and 
various lattice resolutions (Fig. ^ right). Furthermore, the magnitude of the quadratic divergence 
turns out to be C2Ag^Q ^ [30 MeV]^ and essentially the same number follows from the consider- 
ation of the string profile height. Although the above treatment seems to be self-consistent, it is 
crucial to investigate the systematic biases of the ground state extraction, to which we turn next. 



3. Wilson loops 

It seems that the only reliable way to estimate the precision of f — > oo extrapolation is to get 
rid of it and to ensure that only the ground state propagates in the considered channel. To this end 
we combined the multilevel (actually, 3-level) scheme of Wilson loops calculation with the special 



ground state selection technique conventionally known as multichannel QIQQ. In essence it reduces 
to choosing the set {| / )}, / = 1, ...jA^^f of string states and solving the generaUzed eigenvalue prob- 
lem (GEVP) 

Wij {t , R) ef = X W^-y (fo , R) ef . (3.1) 

Here Wij{t,R) is t x R Wilson loop correlator evaluated for the states | / ), \j), which we prepared 
using different number of APE smearing sweeps with the same smearing parameter (see Table[l|). 
The best available approximation to the ground state corresponds to the eigenvector e^"^ with max- 
imal eigenvalue A . For the string profile studies this approach applies literally and leads to 

As={s),-ef\Ws]ijef^/l<^''\ (3.2) 

where [Ws]ij correlation matrix is similar to Wij, but with local action density operator inserted at 
time-like [transverse] distances T [r] from the loop geometrical center (longitudinal shift was kept 
vanishing). Technically, our approach is similar to that of Ref. [^]: at the lowest (deconfining) 
level of our 3-level scheme only the products of time-like open holonomies (possibly with action 
density insertions) were averaged. At the next (confining) level these were accumulated and open 
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E 




R, fm R, fm 

Figure 4: Left: Squared string width 5^ versus distance R plotted in logarithmic scale. Straight line illus- 
trates the bosonic string theory prediction (1.3), C, = 1, with fitted parameter Rq k, 0.2 fm. Right: height of 
the Gaussian profile, As[r^ = 0), as a function of R, sohd hnes represent the rough upper/lower estimates of 
the gluon condensate. 



staple-like segments of Wilson loops were calculated as well (to reduce the statistical errors the 
space-like pieces of the loops were required to be at least /cr/2 distance away from the slab bound- 
aries, l„ being the inverse temperature of deconfinement transition at fixed boundaries). Finally, 
at the upper-most level the averaged open holonomies and their tensor products were combined to 
obtain the correlation matrices Wtj, \Ws\ij. Note that with definition (3.2) As depends upon five 
parameters, {to,t,z,R,r^), where the triple {to,t,T) is used exclusively to ensure the proper ground 
state selection. 

The multichannel approach with three trial states allows to significantly increase the ground 
state overlap. In all cases it was at least 0.99 differing from unity in the third digit well within the 
statistical errors. As a consequence the dependence As{to,t, t) tuned to be trivial as is exemplified 
for j8 = 2.680, 40"^ data on Fig. | (left), where AsatR = 0.93 fm is plotted as a function of T with 
two choices of GEVP window (fo,0 = {(8> 18)> (8i22)} and various transverse distances. The 
straightforward constant fit gives then the t = oo asymptotic of the string profile, which is presented 
on Fig. m, right. To our surprise, the points obtained at various lattice resolutions fall, in fact, on the 
top of each other with no sign of scaling violations. The only conclusion we could draw is that the 
undertaken t oo extrapolations are in fact ambiguous despite of their apparent self-consistency. 
Therefore the previous claim of the string shrinkage is void and at least within the precision of our 
data the quadratic divergence is not seen in the action density difference (1.1). 

Once the claim of linear vanishing of the string width with the lattice spacing had been aban- 
doned it is worth to consider the dependence of 5^ upon the distance R (Fig. ^, left). According 
to (1.3) it is expected to be logarithmically increasing function and indeed Fig. ^ qualitatively con- 
firms this. It is also apparent that the data does not allow quantitative conclusions, only order of 
magnitude estimates could be given. To this end we tried the logarithmic fit of the type (1.3). Al- 
though the data prefer value slightly above one, the difference is insignificant. Thus we fixed 
to unity and evaluated the Rq parameter, the optimal value of which turned to be Rq 0.2 fm. 

Finally, let us remark that for w 1 fm long confining string the height of its action density 
profile is about 0.01 8 GeV"* (Fig. ^, right). Numerically, this is very similar to the magnitude of 
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the gluon condensate in SU(2) gluodynamics [^]. We qualitatively conclude therefore that gluon 
condensate vanishes on the symmetry axis of well developed confining string. 

4. Conclusions 

In this paper we critically reviewed the observation of the confining string shrinkage in the 
continuum limit of pure glue SU(2) lattice Yang-Mills theory. The aim was to address the problem 
of systematic biases coming from the imperfect ground state determination and from the string 
profile distortions due to the presence of singular point-like sources at its end-points. The second 
question is proved to be inessential: in the case of torelon correlator, where the flux tube exists 
without point-like sources, the string shrinkage had been also identified. However, addressing the 
issue of excited states contamination we discovered that it is a very delicate and perhaps not a well 
defined problem of taking the limit of infinite Euclidean times. Despite of apparent 
self-consistency of considered t extrapolations they lead to incorrect results. Namely, in the 
case of Wilson loops it is possible to construct almost perfect approximation to the string ground 
state, the shape of which is in fact independent upon the Euclidean time extent. Then the effect of 
string shrinkage completely disappears proving that it is an artifact of improper ground state 
selection. As a byproduct of our analyzes we qualitatively confirm the validity of the effective 
bosonic string theory description of the string profile and argue that the gluon condensate vanishes 
within the well developed confining flux tube. 
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